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ATTITUDE MOTION OF A NON-ATTITUDE-CONTROLLED 
CYLINDRICAL SATELLITE * 

C. K. WILKINSON | 


ABSTRACT 

In 1985, two non-attitude-controlled satellites were each placed in a low earth 
orbit by the Scout Launch Vehicle. The satellites were cylindrical in shape and 
contained reservoirs of hydrazine fuel. Three-axis magnetometer measurements, 
telemetered real time, were used to derive the attitude motion of each satellite. 
Algorithms are generated to deduce possible orientations (and magnitudes) of each 
vehicle's angular momentum for each telemetry contact. To resolve ambiguities 
at each contact, a force model was derived to simulate the significant long-term 
effects of magnetic, gravity gradient, and aerodynamic torques on the angular 
momentum of the vehicles. The histories of the orientation and magnitude of 
the angular momentum are illustrated. 


• This work was performed in conjunction with Contract F04701-78-C-0125 for the United States 
Air Force , Space Division 

t Principal Engineer , Flight Dynamics Department , Textron Defense Systems 


412 



INTRODUCTION 


The Scout Launch Vehicle placed the two satellites into a low earth orbit inclined 
at 37° to the equator. Approximately one minute after orbit insertion, the 
Scout-4th-Stage/Satellite-System was despun to about 90 deg/s and the two satellites 
were separated sequentially. The separation event imparted lateral rates to the 
vehicles which modified their angular momentums and coning angles The coning 
angles just after the separation event were predicted to be approximately 10° and 
50° for the respective satellites. Each satellite contained a reservoir of hydrazine 
fuel amounting to 15% of the total mass. It was expected that energy dissipation 
from the sloshing fuel would cause the satellites to quickly attain a 90° coning 
angle, i.e. a flat spin about the maximum inertia axis. The satellites are essentially 
axisymmetric with the minimum inertia axis being the axis of symmetry (see 
Fig. 1). The maximum inertia axis should be located very close to the Y c - Z c 
plane. 

PROBLEM STATEMENT 

The problem is to derive each satellite's motion characteristics which are needed 
to validate and interpret satellite system and mission performance. 

ATTITUDE RELATED DATA 

Three-axis magnetometer data are telemetered real time, during contacts, at the 
rate of 8 samples/second. Fig. 2 illustrates the magnetometer data histories (MAGX, 
MAGY and MAGZ along the X , Y and Z axes) for a 7 minute contact with vehicle 2 
on revolution (Rev) 20.4, approximately 1.3 days after launch. The dropouts and 
wild points have not been removed from the illustrated data. The equal periodicity 
on each axis, and the small amplitude in the MAGZ data, indicate that the satellite 
is indeed in a flat spin about an axis very near the Z-axis of the satellite. Both 
vehicles achieved a flat spin prior to the first telemetry contact which occurred 
approximately one hour after launch. 

The telemetry data is a) processed real time for CRT display and associated hardcopy 
output, and b) stored on an analog tape for optional post-pass processing. Both 
processing methods convert the raw magnetometer data to engineering units using 
calibration curves derived on the ground prior to launch. The quantization 
interval, in engineering units, is 0.472 /juT. Real-time hardcopy output of the 
magnetometer data at a rate of 1 sample/second is routinely available. Since post- 
pass analog tape processing is expensive, a ground-rule was established that the 
analysis process be compatible with use of the real-time hardcopy output. 
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ANGULAR MOMENTUM, SINGLE CONTACT 


External torques acting on the orbiting satellites are small. Thus, during a short 
contact, the angular momentum H (referenced to an inertial Newtonian frame) 
may be assumed to be constant. At a given time t, let B(t) be the magnetic 
induction vector, and establish a coordinate system <Xj (0> a 2 (0> and a 3 of mutually 
perpendicular unit vectors satisfying: 
a 3 is parallel to H, 
a 2 (0 is parallel to a 3 x B(t), and 
a 1 (t) = a z (t) x a 3 . 

The geometry at a given instant in time is illustrated in Fig. 3. Let 9 - ( X Y Z ) 

be the principal axes coordinate system of the vehicle with Z the maximum 

inertia axis and X the minimum inertia axis. For flat spin motion, Z is parallel 
to H, i.e. Z = a 3 . The relationship between the vehicle magnetometer coordinate 
system Ji = (X M f M Z M ) and 9 is illustrated in Fig. 4. The magnetometer axes are 
designed to be mutually orthogonal and oriented along the vehicle geometric, 
frame. The magnetometer data gives no information on the frame. We thus 
work with the M frame and assume it is an orthogonal system. 

The vehicle motion in principal axes is described by: 

X(l) = &i(T 0 ) cos (Q p t + i/J + a z (T 0 ) sin (0 p t + ^ 0 ) 

f(0 = -0,(7;) sin (cy + ^ 0 ) + a 2 (7;)cos(cy + ^ 0 ) d) 

Z(t) = a 3 

where 

O p is the precession rate, 

^ is the initial precession angle, 

T 0 is the initial time, 

T is the current time, and 

i = T - T 0 . 

In the M frame, making the small angle approximation for 6 and setting cos 5 
to 1: 

X M (t) = X(t) + <5 sin 0, f(0 - <5 cos 0^ Z (0 

Y M (t) = 6 sin <f> 3 X(t) + cos $ f(0 + sin C Z(0 (2) 

Z M (t) = <5cos 0 3 X(O - sin^T(0 + cosfZCO 

where 

0 6 and 0 3 are the Euler angle rotations about the X, Y and X axes, 

respectively, and 
? = 0j + 0 3 (See Fig. 4.) 
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The vehicle's spin period is short relative to its translational motion. Thus we 
— ► 

assume that B, and therefore a x and a 2 , are constant over one spin cycle. Then 
the magnetometer history over one cycle satisfies: 

X M B = \b\ [- <5 cos At cos 0j - 6 sin At sin 0 X sin(O p / + '0 o ) 4 sin At cos{Q p t 4 yj] 

? M B = \B |[ cos a* sin f - sin At cosCsin(Q p < 4 yp Q ) 4 <5 sin /x sin 0 3 cos(Q p < 4 yj] (3) 

4 B = |i?| [cos /t cos f 4 sin /x sin f sin(O p < + -0 O ) 4 <5 sin ix cos 0 3 cos(O p / 4 -0 O )] 

The extrema, fT 1 , in (3) are. 


II 

|i?| [ ± sin ix - 

6 COS IX COS 0 X ] 

(4a) 

II 

V £ 

|f?| [ ± sin [x cos £ 4 

cos At sin f ] 

(4b) 

II 

|i?| [ ± sin At sin f 4 

cos /z cos f ] 

(4c) 

The extrema, E ± , in the 

observed values, assuming no scaling errors, are: 

£* = £* +3 

( 5 ) 


where 


D = 


D* 

By 


are the biases along the magnetometer axes. 


14 #/ 

The derivatives of (3) are: 


dX M B 
di 

dY M B 

dt 

d Z M ■ B 

dt 


= - 0 p | f? | sin /J. [ 6 sin 0j cos(O p / 4 -0 o ) 4 sin(O p / 4 t/' 0 ) ] 


= - 0 p |i?| sin a* [ cosf cos(O p / 4 -0J 4 5 sin 0 3 sin(O p / 4 -0 o )] 


Op |/?| sin fx [ sin f cos(O p / 4 Tp Q ) - 6 cos 0 3 sin(O p / 4 i p o ) ] 


The occurrences of the extrema are summarized below: 


( 6 ) 


4 

• B 

: tan(O p / 4 -0 O ) = 

- <5 sin 0 X 

(7a) 

4 

B : 

tan(O p / 4 -0 O ) = 

- cos £ / <5 sin 0 3 

(7b) 

4 

B : 

tan(O p / 4 -0 O ) = 

sin f / <5 cos 0 

(7c) 


Equation 7c is particularly useful in determining 6 when £ and 0 3 are small. In 
this case C w 0 r Note that Equation 7a is w 0. 
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Procedure 


1) Obtain /z from the amplitude in 4 B (Equation 4a). 

2) Obtain f from either 4 B or 4 B , whichever has the smaller 
amplitude (Equation 4b or 4c). 

3) Equations 7b and 7c yield approximations to <5 and 0 3 from inexact 
values of t determined from the magnetometer data. 

4) Determine E ± from Equation 4. 

5) Determine D from Equation 5. 

If sin ^ 1, then an accurate value of /x or D cannot be obtained. The best 

procedure here is to assume sin /x = 1, perform steps 2) and 3), and use an average 
— * 

value of D from other contacts to solve for cos/z from either of Equations 4b 
or 4c. With the exception of this case, the magnetometer biases need not be a 
priori known to solve for /x and 

Determination of Quadrant 


The following rules resolve the quadrants of f and /x. 


sgn ( sin f ) 

sgn ( cos <■ ) 
sgn ( cos ix ) 


( d X M B\ 

sgn I — I at t maximizi 

/ d X M B\ 

sgn I -77 j at t maximizi 

sgn ( Z M B ) sgn ( sin f ) . 


y m b. 

Z M B . 


( 8 ) 


Effect of B Variation with Time 


In practice, fixing B has little effect on the solutions for ix, f and <5. The variation 
in Bit) is important in the determination of Q p : 


Q 


v 


°p + 


d aXt) 
dt 


0 ) 


Specifically, the observed precession rate, Q p , can be quite different from 0 p if 
—► — ► — ► 

H ■ B is small. Further, the values of a 2 (t ) for the two solutions of H are opposite 

in sign. 
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Example 1: sin /x is Small 


Table 1 gives results for the contact on Rev 164 with vehicle 2. The telemetered 
magnetometer histories are illustrated in Fig. 5. The B values were obtained from 
an orbit simulation code using a 12 degree spherical harmonic expansion for the 
geomagnetic potential. The algorithms are used to compute ytx, (£ + £ )/2, D 
and <5. Independent estimates, D, for the biases, for comparison purposes, were 
obtained by an estimation scheme which minimizes, in a least squares sense, the 
difference between the magnitudes of the bias-adjusted measured /vector and the 
modelled magnetic induction vector. The value of the biases vary with the 
operating configuration of the satellite which is slightly different for each of the 
three times The accuracy in <5 is poor because of the small variation in Z M B 
and the large quantization interval. The offset in the extrema of Z M B from 
B and y M ■ B is illustrated in Fig. 6 for Rev 20.4 which gives a larger variation 
in ■ B The satellite rotates about 64° between the occurrences of the X M • B 
and Z M ■ B maxima. This angle was used in Equation 7c to solve for 6 assuming 
0 3 is 0. Analysis of several contacts indicated that the values for £ and 6 were 
close to 5° and 2°, respectively. When used in Equations 4 and 5, these yield a 
value of D which agrees very closely with D 


Example 2: sin fji w 1 

It is difficult to determine /x accurately when sin n w 1. Table 2, pertaining to 
the contact on Rev 20.4 of vehicle 2, illustrates the procedure described above. 
(Refer to Fig. 2 for the magnetometer histories.) First obtain estimates, jx and D, 
of ytx and D Equation 4a is used to determine /x. In general, D is an average over 
several selected contacts. (The value of D used in Table 2 was obtained using the 
aforementioned bias estimation code since the analog tape for that contact was 
processed post- pass.) Then E i is determined from Equation 5 and ix is computed 
more accurately from (4c) in the form: 


cos 


£+ 4 £~ 

2 | £ | cos C 


where £ is obtained by the usual procedure, but using The other components 
of D now agree fairly well with those of D. 
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TABLE 1 COMPUTATIONS FOR REV 164.2, VEHICLE 2 



T = 

82849s, 

GMT 

T = 83029s, 

GMT 

T = 

83179s, GMT 

X 

Y 

Z 

X 

Y 

Z 

X 

Y Z 


-5.2 

10.6 

38.2 

-10.3 

7.4 

43.0 

-6 6 

11.3 45 3 

E~ 

-22.8 

-7.0 

36.3 

-17.0 

0.7 

42.3 

-20.6 

-2.6 44.0 

131 


27 62 



30 39 



32 92 

C 


6.2 



6.0 



5 3 

6 


3.0 



2.9 



2.6 



18.6 



6.3 



12.3 


-1.36 

2.82 

26.02 

-1 52 

3.16 

30.04 

-1.46 

2.97 32.03 

D 

-12 64 

-1 02 

11 23 

-12.13 

0.89 

12.61 

-12.15 

1.38 12 62 

D 

-13 17 

0.15 

11.30 

-12.62 

1.61 

12 55 

-12.62 

1 61 12.55 


TABLE 2 

COMPUTATIONS FOR REV 20.4, 

VEHICLE 2 


T = 

37941s, 

GMT 

T = 38147s, 

GMT 

T = 

38327s, GMT 


X 

Y 

Z 

X 

Y 

Z 

X 

Y Z 

E+ 

20.4 

33.7 

16.1 

20.3 

34.6 

8.5 

17.0 

30.7 0.8 

E~ 

-47.4 

-33.3 

9 8 

-45.7 - 

-30.8 

2.3 

-42.7 

-28 4 -48 

I3| 


34 09 



33 92 



33 32 

D 

-13.53 

0.17 

12 48 

-13.08 

2.64 

13.47 

-13.20 

2.32 13 05 

A 


90. 


103.4 



116.4 


0.03 

0.03 

0.47 

0.38 

-0.74 

-8.07 

0.45 

i -1.17 -15.05 

? 


5.3 



5.4 



5.3 

6 


2.6 



2.6 



2.6 



89.2 


103.8 



117.0 

— ► 
D 

-13.47 

0.15 

12.48 

-13.09 

2.66 

13.47 

-13.43 

2.55 -15.05 
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Fig. 5 Magnetometer Histories for Rev 164.2, Vehicle 2 



Fig. 6 Magnetometer Cycle Offset for Rev 20.4, Vehicle 2 
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Ambiguity in H 


— Y — ► 

- The procedure described above resolves H to lie on a cone of angle /x about B . As 

the geometry, specifically B, varies with time due to satellite orbital motion, H 

can be resolved to two points: the intersection of two cones. Over a short contact, 

the magnetic induction vectors are essentially co-planar. The two possible H 

vectors are the true H and its mirror image relative to the 'B ( t ) plane'. The key 

point is that H cannot be uniquely determined over a short contact. This can be 

observed very readily in Fig. 7 which illustrates the loci of possible H (minor 

— ► 

circles of radius /x i about point Bi which is the intersection of B i with the unit 
sphere) in right ascension and declination at the three times in Table 1. (A Kalman 
filter estimation technique designed specifically for this problem could not, in 
general, decide which of the two possibilities was the correct H orientation.) 

The geometry of the two possible H solutions will vary with the orientation of 

the actual H with respect to the B vectors during the contact. For Rev 21.4, 

vehicle 1 (Fig. 8), H is quite far from the 'plane' of the B(t ) vectors so that the 
— ► 

two H possibilities are widely separated. This contact is sufficiently long - 7 minutes 
- for B(t) to be non-planar which allows the correct H to be resolved. On Rev 20.1, 
vehicle 2 (Fig. 9), H is many degrees from B but is close to the B(t) plane so that, 
with the uncertainties associated with /n, the possible H values lie along an arc 
of « 30° encompassing declinations from -20° to -45°. Rev 1.1 for vehicle 1 (see 
Fig. 10) has similar geometry. 

Since the satellite contacts are short, on the order of 3 to 5 minutes, and may 
^ be separated by about 30 days, other methods must be employed to resolve the 
ambiguity in the angular momentum orientation. The next section discusses how 
this issue is resolved. 

ANGULAR MOMENTUM RESOLUTION 

Resolution of the correct angular momentum orientation between the two 
possibilities requires a) an initial orientation for H , and b) a torque model from 
which a long-term history of H can be derived. 

Initial Orientation of H 

There were frequent contacts with the satellites during the first few days in 
orbit. Fig. 10 illustrates a resolved orientation for the H of vehicle 1 at the 
intersection of the loci for contacts on Revs 1 and 2. This solution is essentially 
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on the mean of the expected loci determined by a statistical analysis, conducted 
pre-launch, of the separation of the satellites from the 4th Stage of the launch 
vehicle. Fig. 10 illustrates a likely H path from this solution to the H solution 
on Rev 21 This corresponds to a precession in H of « 15° /day 

The initial orientation of H for vehicle 2 was established at a right ascension of 
162° and a declination of -30° 

Torque Model 

Torques acting on the orbiting satellites originate from magnetic, gravity gradient, 
aerodynamic, and solar radiation pressure effects. The approach taken was to 
obtain expressions, available in the published literature, for the effect, on H , of 
each significant torque over one satellite revolution. These expressions, along with 
a dipole model of the earth's magnetic field 5 were incorporated into a semi-anlytic 
orbit generation code using a 1-rev step-size. The torque expressions are described 
briefly below. Because of the small size of the satellites, solar radiation pressure 
was judged to be a relatively insignificant contributor and was not included in 
the model. 

Magnetic Torques: - Magnetic torques result from the interaction between the 
magnetic properties of a spacecraft and the ambient magnetic field of the earth. 
The primary magnetic disturbance torques are: 

1. Dipole moment from the permanent magnetism in the spacecraft; 

2. Eddy currents induced when a conducting body moves in a magnetic 
field; 

3. Spacecraft generated current loops; and 

4. Hysteresis damping. 

Because of the limited amount of spacecraft operating time and the specific 
amount and properties of permeable material present on the spacecraft, items 3. 
and 4. were judged to have small contributions to total vehicle torque and were 
not included in the model. 

The torque, T D , due to the dipole moment is normal to a 3 
has only a precession component. It satisfies 2,3,5 : 

T D = m d a 3 x B 

where 

m d is the satellite's dipole moment along its Z-axis. 


= h/\ H I, and hence 

(10) 
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Making the assumption that the satellite Z-axis is inertially fixed for an orbital 
period, the average torque over one orbit rev satisfies: 

(T d ) av = m d“3 x (5) av (11) 

(i?) av can be integrated with respect to time over an orbital period 2 to give an 
average induction vector. 

The total torque due to eddy current effects satisfies 5 : 

T E = k e (u z xB)xB (12) 

where 

k e is a constant which depends on the geometry and conductivity 
of the rotating object, and 
co z is the angular velocity vector of the satellite. 

Equation (12) can be separated into despin, Tez' and precession, T ev components 
of torque. 

T ez = - k e {B L fu z (13) 

T'eI = ^e^z^l^z ( 14 ) 

where 

B , is the component of B orthogonal to Z, 

B z is the component of B parallel to Z, and 

c o z is the angular velocity. 

These equations are conceptually simple. The complication arises in finding the 
average values for {B^f and B z B i over one revolution of the satellite. Reference 3 
derived equations for B^ and B z and then integrated the resulting products. It is 
easier to a) resolve (« z x b} x B into its three components, b) integrate, and then 

c) compute the despin and orthogonal components, (T EZ )^ and (7’ ZI ) av . of torque. 

Since the orbit eccentricity is small, terms of o(e 3 ) and o(e 4 ) can be deleted. 


Gravity Gradient Torque: - The average torque over one orbit due to the effect 
of the earth's gravitation on the satellite is 1,2,3,4 : 



1.5 - 3 — - V n i ( 
a 3 (l - e 2 ) 1 5 ' 


^i)(a 3 N)(a 3 x N) 


where 

is the earth’s gravitational constant, 
a is the semi-major axis of the orbit, 

e is the eccentricity of the orbit, and 

N is the unit normal to the orbit plane. 


(15) 
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Aerodynamic Torque: - The average aerodynamic torque over one satellite spin 
cycle satisfies: 


T A = ±pv z SC D d(l v X d 3 ) (16) 

where 

p is the density of the atmosphere, 

v is the velocity of the satellite, 

S' is the projected area in the direction of motion averaged over 
one satellite rotation, 

C D is the drag coefficient, 

d is the distance from the center of mass of the satellite to its 
aerodynamic center of pressure along Z averaged over one satellite 
rotation, and 

I v is the unit velocity vector of the satellite. 

The average torque over one orbit was determined using a relatively common 
procedure 3 for somewhat elliptical orbits. 

Spin-axis Rate Change in One Orbit Revolution: - Assuming that the vehicle is in 
a pure spin motion about the body Z-axis, i.e. the principal inertia axis, the 
change in the spin rate in one orbit revolution caused by the eddy current torque 
is: 

~ ~ | (Trz) a v | 

where 

I z is the moment of inertia about the body Z-axis, and 
P is the satellite orbital period. 

Spin-axis Precession in One Orbit Revolution: - Over one orbit rev, the angular 
momentum vector, H, precesses by the angle |r av |/ <o z I z from H towards T Ay 
where 7 av is the sum of the individual torque contributors normal to H . 

RESULTS 

With the incorporation of the torque model into the orbit generation code, the 
determination of approximate values for the key torque-related vehicle-unique 
properties proceeded rapidly using an iterative simulation process. The inertia 
properties given in Fig. 1 and an aerodynamic lever arm, d, of 0.75 cm, equal to 
the pre-flight measured e g. -offset along the Z c -axis, were used initially for both 
vehicles. Though they were varied to determine sensitivities, these values were 
used in the final simulation. The derived magnetic properties are summarized 
in Table 3. 
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TABLE 3 DERIVED MAGNETIC PROPERTIES 


VEHICLE 

DIPOLE MOMENT, m H 



A - m 2 

mVohm 

1 

+0 543 

465 

2 

-0.230 

465 


The dipole moment values are within expectations since satellite magnetic properties 
were minimally controlled during the design process. The k e constant should be 
the same for both vehicles The derived value is reasonably close to the 600 
m 4 /ohm value computed pre-flight. The histories of the orientation of H over 
the first six months in orbit are illustrated in Figures 11 and 12 for vehicles 
1 and 2 respectively. Over this period, the maximum differences between the 
simulated and calculated H orientations were one day ( w 15 revs) in-track (parallel 
to the trace of H), and w 10° cross-track (normal to the trace of H). Generally 
the agreement was much better. Fig 13 illustrates, for Rev 3056 of vehicle 2, 
typical agreement between the simulation and computations for the orientation 
of H. 

The simulated spin rate histories are summarized in Figures 14 and 15. They 
match the 'observed' data very well except for the first few days where the 
differences are attributed to residual angular momentum in the hydrazine fuel. 
There appears to be no significant hysteresis in the material magnetization cycle. 
This would be manifested as a constant component to the primarily exponential 
spin decay rate. 

Figures 16 and 17 illustrate the histories of the angles between H and N . For 
vehicle 1, H is never more than 20° out of the orbit plane. Its projected area 
along the flight path, averaged over a spin cycle, is essentially constant. For 
contrast, the H for vehicle 2 is at times nearly normal to the orbit plane. Thus, 
its average projected area along the flight path is smaller than that for vehicle 1. 
This correlates with observed periods of w 5% smaller effects of drag on vehicle 2. 

Finally, Figures 18 and 19 illustrate the dynamic history of the sun aspect angle. 
This, along with the eclipse history, are necessary inputs to an assessment of the 
thermal performance of the vehicles. The long thermal time constant («*9 days) 
of the well insulated vehicles, combined with the rapid fluctuation in the sun 
aspect angle, led to a very uniform internal temperature history. 
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Fig. 12 Evolution of the Orientation of H for Vehicle 




SUN/SPIN-VECTOR ANGLE, degrees (Q* ORB IT-NORMAL/SPIN- VECTOR ANGLE, degrees PRECESSION RATE, degrees/second 



14 Precession Rate History for Vehicle 1 Fig 15 p reC ession Rate History for Vehicle 2 
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